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Abstract. XMM-Newton observations of 19 galaxy clusters are used to measure the elemental abundances and their spatial 
distributions in the intracluster medium. The sample mainly consists of X-ray bright and relaxed clusters with a cD galaxy. 
Along with detailed Si, S and Fe radial abundance distributions within 300-700 kpc in radius, the O abundances are accurately 
derived in the central region of the clusters. The Fe abundance maxima towards the cluster center, possibly due to the metals 
from the cD galaxy, are spatially resolved. The Si and S abundances also exhibit central increases in general, resulting in 
uniform Fe-Si-S ratios within the cluster. In contrast, the O abundances are in general uniform over the cluster. The mean O 
to Fe ratio within the cluster core is sub-solar, while that of the cluster scale is larger than the solar ratio. These measurements 
indicate that most of the Fe-Si-S and O in the intracluster medium have different origins, presumably in supernovae la and II, 
respectively. The obtained Fe and O mass are also used to discuss the past star formation history in clusters. 
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1. Introduction 

Metals in the universe are being produced as a product of stellar 
evolution and ejected into intergalactic space via stellar winds 
and supernovae. The location of these metals can be mea- 
sured by X-ray observations of the intra-cluster medium (ICM). 
These observations indicate that the ICM contains an amount 
of metals comparable to the total amount of metals found in 
galaxies (e.g. Tsuru fTOOTI Arnaud et al. 1992). Therefore, X- 
ray measurement of the ICM metallicity distribution is a unique 
way to study the chemical history of galaxies. Most impor- 
tantly, the total amount of metals and the abundance ratios 
among elements (e.g. O/Fe ratio) constrain the total number 
of stars that have formed as well as their origin (i.e. the relative 
ratios of supernova types; Renzini et al. ll993> . 

Measurements from ASCA and BeppoSAX observations re- 
vealed important properties of the ICM metallicity. These in- 
clude an excess of Fe around cD galaxies (e.g. Fukazawa et al. 
|2000 > De Grandi and Molendi|200D and variations of the Si/Fe 
ratio within a cluster (e.g. Finoguenov et al. 2000) and among 
clusters (e.g. Fukazawa et al. I1998> . The variations in Si/Fe 
suggest that the metals in the ICM have at least two different 
origins, presumably SN la and SN II. The spatial resolution 
of these previous instruments was limited, causing the derived 
total amount of metals to depend on their assumed spatial dis- 
tribution. The low spatial resolution combined with the limited 
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spectral resolution also introduced a systematic uncertainty in 
the temperature structure. This is particularly important in the 
central cool regions. This ambiguity in turn caused severe er- 
rors in the abundances of Fe and other elements. Furthermore, 
in most cases these measurements are limited to the Fe, Si and 
S abundances. 

Recent observations with Chandra and XMM-Newton have 
allowed a detailed study of the ICM metals, using an improved 
spatial resolution and sensitivity not only for the Si, S, and Fe 
lines but also for the O emission. For example, Tamura et al. 
(2001 bi measured the spatial distribution of the O, Si, S, and 
Fe in Abell 496 and found a difference between O and Si-S- 
Fe. David et al. i20Q\i used the Chandra data of the Hydra- 
A cluster to measure the Si and Fe abundances. In addition, 
the abundance distribution of several elements in the central 
region (70 kpc in radius) of the Virgo cluster was obtained with 
high accuracy (Finoguenov et al. 120021 Gasteldello & Molendi 
I2OO2I Matsushita et al. l2003t . 

In this paper, we use the XMM-Newton EPIC and RGS 
instruments and measure the spatial distribution of heavy el- 
ements in 19 clusters. This systematic analysis enables us to 
compare the metal distribution in the sample. In addition, we 
may reduce the statistical eiTors on any given measurement 
when the variation in the sample is negligible. We selected the 
sample based on their X-ray flux and spatial size. Therefore the 
sample is limited to nearby X-ray bright clusters with a mean 
redshiftof0.06. 
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We obtained the radial distribution of Fe on a sub-arcmin 
scale for each cluster and confirmed the central increase of the 
Fe abundance within a typical radius of 100 kpc in most cases. 
Furthermore, radial profiles of O, Si, and S relative to Fe have 
been measured in a fair number of clusters. We found that the 
O and Si-S-Fe abundances have different spatial distributions 
when we combine the results from several clusters. 

Throughout this paper, we assume the Hubble constant to 
be Hq = lOO/i kms^^ Mpc~^ and use 90% confidence lev- 
els unless stated otherwise. Solar abundances are taken from 
Anders & Grevesse i fT989b with Fe/H=4.68 x IQ-^. 

2. Observations 

We analyzed 19 clusters taken from GT and PV observations 
of XMM-Newton (TableQ. The RGS exposure times are some- 
times shorter than those reported in Peterson et al. (2003), be- 
cause of a difference in the criteria for rejecting the exposure 
period with high background. We have previously reported de- 
tailed results for A 1795 (Tamura et al. 2001ai, Sersic 159-03 
(Kaastra et al. 2001 ), A 1835 (Peterson et al. l200Tl . and A 496 
(Tamura etal. 2001bi. 

The present paper is one of a series of papers with a sys- 
tematic study of cluster properties using similar data sets: soft 
X-ray excess (Kaastra et al. '2003a'); central cooling region us- 
ing RGS data (Peterson et al. 2003 1; central cooling using EPIC 
(Kaastra et al. 2004 1. Here we focus on the spatial distribution 
of the metal abundances using both the EPIC and RGS data. 

A detailed description of the XMM-Newton satellite and 
its instruments is found in Jansen et al. (2001), Turner et al. 
(EPIC/MOS; 2001), Struder et al. (EPIC/pn: l200lJ . and den 
Herder et al. rRGS: .200lJ . 

3. The EPIC analysis 

3.1. Spectral extraction 

In order to derive the spatial metallicity distribution, we used 
spatially-resolved EPIC spectra. We used the Science Analysis 
System (SAS; version 5.3.3) for the basic data processing. 
Then, we corrected the filtered events for the background, the 
telescope point spread function (PSF), and the effective area. 
We described the data correction and extraction extensively in 
Kaastra et al. i2004l . 

We removed high background time periods based on the 
number of counts in the 10-12 keV band over the full field 
of view. Even after rejecting these high background events, 
there were still variations in the instrumental background from 
observation to observation. In order to correct for thiso weak 
variation, we sorted the cluster and blank sky data into sub- 
sets with the same 10-12 keV count rate. Then, the frequency 
distributions of source count rate and a set of the blank sky 
spectra at different 10-12 keV count rate level were used to de- 
termine the subtracted background spectrum. The source and 
background spectra were obtained from the same detector re- 
gions. To compensate for remaining variations in the instru- 
mental background, we included a systematic uncertainty of at 
least 10 % in the spectral fitting. 



In both MOS and pn the Al lines around 1 .48 ke V, which 
are close to the line positions of H-like and He-like Mg ions, 
are strong. Therefore, we did not attempt to derive the Mg 
abundance except for bright cluster central regions. In addition, 
in the MOS, there are instrumental Si lines around 1.74 keV. 
Compared to the Al lines, these lines are relatively weak and 
not as close to that from astrophysical Si ions. Therefore, an 
additional uncertainty in the Si background lines produces no 
change in our result. We verified this using representative clus- 
ter data. 

In the low energy band, the cosmic background dominates 
the instrumental background and varies spatially by more than 
10%. Therefore, we included larger systematic errors in the soft 
X-ray band as shown in Tabled 

In all observations the bright X-ray point sources were re- 
moved. Typically, a median of 8 point sources (for MOS 1 and 
MOS2) or 15 (forpn) were discarded. 

Then, we made a background subtracted radial profile for 
each energy. This profile was corrected for the exposure, tele- 
scope PSF, and vignetting. We assumed energy-independent 
PSF correction (Kaastra et al. I2004> . Therefore it affects only 
the spectral normalization, not its shape. On the other hand, the 
correction for the vignetting due to the telescope is dependent 
on the X-ray position and energy. The MOS cameras also have 
a correction due to obscuration of the Reflection Grating Array 
that is dependent on position and energy as well. 

From the corrected radial profile for each energy band, we 
extracted projected spectra from concentric annuli centered on 
the X-ray maximum with outer radii of 0.5, 1, 2, 3, 4, 6, 9, 12, 
16 arcminutes. In addition, assuming spherical symmetry, we 
derived deprojected spectra from shells with the same radial 
binning as above. The deprojected spectrum represents a set of 
count rates in a spherical shell between the two spheres and 
can be given by a linear combination of the observed number 
of counts in the (projected) surface brightness profile. Note that 
we ignored X-ray emission outside the field of view and that the 
spectra of the outermost annulus (12'-16') are used only for the 
deprojection for the inner shells. We used different innermost 
and outermost bins among our clusters as shown in Table ^ 
depending upon the statistical quality of the spectra. 

We do not completely take into account the spectral mixing 
due to the PSF. Therefore the true profiles of temperature and 
metallicity could be much steeper than the present results. This 
effect should be small since the the radial bin size is larger than 
the PSF width (a half energy width of 15"). 

3.2. Spectral fitting 

In most cases, the instrumental background dominates the 
source emission above ^ 8 keV. Between 7-9 keV, the 
pn instrumental background contains strong spectral lines. 
Furthermore, when we fitted the pn spectra of bright clus- 
ters with representative models, we found systematic residu- 
als below the Fe-L lines. This is probably due to uncertain- 
ties in the energy redistribution response at the low energy. 
Taking this into account, we restricted our analysis for the MOS 
and pn spectra between 0.3 and 8 keV and between 0.7 and 
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Table 1. Cluster sample. 



Cluster redshifF kT^ r^iT^ C^" XMM Exposure time (ks)'= ¥ 







(keV) 


(10^** m-2) 




revolution 


MOS 


pn 


RGS 


IN 


OUT 


NGC 533 


0.0175 


1.3 


3.00 


C 


195 


38 


31 


49 





6 


A 262 


0.0155 


2.2 


8.94 (5.33) 


C 


203 


24 


15 


35 





12 


Sersic 159-03 


0.0572 


2.4 


1.79 


C/SE 


77 


31 


28 


36 





9 


MKW9 


0.0402 


2.6 


4.18 


C 


311 


22 


13 







3 


2A 0335+096 


0.0344 


3.0 


28.71 (18.64) 


M 


215 


11 


5 


18 





9 


A 2052 


0.0356 


3.1 


2.91 


M/SE 


128 


30 


23 


28 





12 


Hydra A (A 780) 


0.0550 


3.4 


4.80 


M 


183 


17 


3 


25 





6 


MKW3s 


0.0455 


3.5 


2.89 


M/SE 


129 


37 


32 


_ 





12 


A 4059 


0.0466 


4.0 


1.06 


M 


176 


34 


6 


44 





9 


A 1837 


0.0707 


4.4 


4.38 


M 


200 


48 


40 


„ 





6 


A 496 


0.0322 


4.4 


6.44 (4.23) 


M 


211 


9 


5 


18 





12 


A3112 


0.0756 


4.5 


2.61 


M/SE 


191 


22 


16 


24 





9 


A 1795 


0.0639 


5.8 


1.01 


M/SE 


100 


29 


22 


37 





9 


A 399 


0.0706 


6.2 


10.90 


H 


127 


13 


6 




0.5 


9 


Perseus (A 426) 


0.0179 


6.5 


14.90 


H 


210 


43 


33 







12 


A 1835 


0.2541 


7.2 


2.32 


H 


101 





25 


29 





6 


Coma (A 1656) 


0.0240 


7.5 


0.89 


H/SE 


86 


15 


13 




0.5 


12 


A 754 


0.0561 


8.0 


5.67 (4.37) 


H 


262 


14 


11 




0.5 


12 


A 3266 


0.0614 


8.7 


1.60 


H 


153 


24 


19 




0.5 


12 


Virgo^ 


0.0027 


2.6 


1.80 


C 


97 


34 


25 
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^ Redshift, to be used for the distance estimate. 

The temperature outside the central cool region obtained from our EPIC analysis (Kaastra et al. '2004V 

Galactic column density to be used for the spectral modeling. In most cases, these are obtained from H I map (Dickey & Lockman ll990> . 
In other cases, we present the used (X-ray) value with the H I value in parentheses. 

Classification of clusters used in the present paper. "C", "M", and "H" indicate cool, medium-temperature, and hot clusters, respectively. 
"SE" indicates clusters which show soft X-ray excess. 

Net exposure time, after rejection of high background time periods. We did not use the RGS data for clusters indicated by because of 
poor statistics. 

^ Innermost and outermost radii for the EPIC spectral analysis in arcmin. 
^ The data for the Virgo cluster are used only for the EPIC calibration. 



7.0 keV energy ranges, respectively. The three spectra (MOSl, 
MOS2, and pn) were fitted simultaneously with the same spec- 
tral model, but we allowed the relative normalizations to be free 
parameters. 



Table 2. Systematic errors relative to the flux. These were used 
for the EPIC spectral modeling and error estimates. 



Energy band (keV) 


0.3-0.5 


0.5-0.7 


0.7-2.0 


2.0-8.0 


source 


5 


5 


5 


10 


background 


35 


25 


15 


10 



We used a collisional ionization equilibrium model (CIE; 
implemented in the SPEX package' ;Meweetal. 1985 Kaastra 
119921 Liedahl et al. IT995t modified by photoelectric absorp- 
tion. The abundance ratios of C and N were fixed to the solar 
values, based on the best-fit value derived from the RGS in the 
central region of M87 (Sakefliou et al. l2002l l. Those of O, Ne, 
Mg, Si, S, Ar, Ca, and Fe were left free, while those of other 
heavy elements were fixed to that of Fe (e.g., Ni/Fe=l). Thus 

' see ,http://www.sron.nl/divisions/hea/spex/| 



the free parameters are the normalization, temperature (kT), 
and metal abundances (O, Ne, Mg, Si, S, Ar, Ca, and Fe). 

We used a fixed column density for each cluster as shown 
in Tabled These column densities have been adjusted to match 
the observed EPIC spectra as follows. In most cases, the best- 
fit value is consistent with the Galactic one obtained from 
the H I map (Dickey & Lockman II 9901 1 . In those cases we 
kept the column density fixed to the Galactic value. In A 262, 
2A 0335H-096, A 496 and A 754, however, the data require 
additional absorption, possibly caused by Galactic gas and 
dust (we investigated these individual cases in more detail in 
Kaastra et al. ,2003a i. 

Usually, the best-fit parameters (normaUzation, kT, and 
abundances) are derived at the same time by fitting the spec- 
trum in a wide energy band. If however the response matrix 
is not accurate due to a calibration error in some wide energy 
band of the spectrum, the modeled local continuum may dif- 
fer slightly from the actually observed continuum in that wide 
energy band. This could result in errors on the derived equiv- 
alent width of lines in that band and hence in errors on the 
obtained abundances. In fact, when we used this global fitting 
method with the MOS and pn spectra separately with bright 
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clusters, we found that the MOS gives systematically higher 
Si and S abundances by ~ 0.2 of the solar value than the pn 
does. To avoid such errors as much as possible, we derive the 
parameters in two steps. First, we measure the best-fit kT by 
fitting the wide energy band spectra. Secondly, we fix kT, use 
the local spectra in the limited energy band close to the rele- 
vant elemental-ion lines and derive the abundance by allowing 
the normalization to be free. We used energy bands of 0.35- 
1.0 (for O), 0.33-2.2 (Ne and Mg), 1.2-3.2 (Si and S), 2.7^.7 
(Ar and Ca), all in keV. An exception is the Fe abundance, for 
which we use the wide energy band to include both the Fe-L 
and Fe-K lines simultaneously. 

We analyzed both projected and deprojected spectra. The 
projected spectra in general have better statistics but poorer fits 
for each radius, mainly because of the contamination from dif- 
ferent radii with different kT. Also, the derived spatial distribu- 
tion of the temperature changes by using deprojection, in par- 
ticular in the inner regions (See Fig. 2 in Tamura et al. 2001b 
for an example). To obtain an accurate abundance distribution, 
proper modeling of the temperature distribution is necessary. 
Therefore, we present the results of the deprojected spectra. 

3.3. Systematic uncertainties 

We used several calibration sources to estimate the uncertainty 
in the response calibration. Based on this, we place conserva- 
tively systematic errors on the source spectra. Table |2] sum- 
marizes the values of these systematic errors that we used. 
These errors along with those of the background were added 
in quadrature in the binned spectra. 

To assess additional systematic errors, we derived abun- 
dances (Si, S, and Fe) separately from the MOS and pn spectra 
in representative clusters (M 87, A 262, and A 1795). For the 
Fe abundance determination, we separately used Fe-L and Fe- 
K lines. The results are shown in Fig. [2 The pn detector still 
gave systematically smaller Si and Fe abundances than MOS. 
The differences were less than 10-20 % of the obtained val- 
ues and less than the typical statistical errors for each spec- 
trum. In other cases, the two instruments provided consistent 
results. Because the Ovill Lya line was contained only in 
the MOS spectrum in our case, we cannot check the O abun- 
dance in a similar way. Nevertheless, we confirmed that the O 
abundances obtained from MOS are consistent with those ob- 
tained from RGS as shown below (Fig.|2j. The three types of 
XMM-Newton instruments have different energy responses and 
have been calibrated independently, except for the telescope 
response. Therefore the above consistency checks indicate that 
there is negligible additional systematic error above the values 
that we already assigned compared to the statistical uncertainty 
in our sample. 

3.4. Results 

In general, the single temperature (IT) model provided a good 
fit to the spectra. Figures |3] and 0] show examples of the fit. 
However, the spectra from the following regions cannot be de- 
scribed by the IT model (x^/v > 1.2). NGC 533 (0' - I'.O), 



ID 



a 



- 0/ 


Fe 
































































A1B35 







0.5 1 1.5 

From EPIC (r<l') 



Fe/H 




aA0335 



0.5 1 

From EPIC (r<l') 

Fig. 2. Comparison between the EPIC and RGS results for O/Fe 
and Fe/H. 



A 262 (0' - 0'.5), 2A 0335H-096 (C - I'.O), MKW 3s (6' - 12'), 
A 2052 (0' - 0'.5,9' - 12'), A 4059 (C - 0'.5), A 3112 
(1' - 2', 3' - 6'), and Perseus (C - 2', 3' - 4'). 

In the inner regions of some clusters (r <1'), the data ex- 
hibit a clear deviation from the IT model. A major part of this 
discrepancy can be removed by adding an additional cool com- 
ponent. In addition, some spectra in the outer regions cannot 
be modeled by the IT model, probably because of soft excess 
emission below 0.5 keV. We reported earlier a detailed anal- 
ysis of this soft excess in Sersic 159-03, A 2052, A 1795, 
and Coma (Kaastra et al. I2003at . In addition to these clus- 
ters, A 3112 also shows a similar excess above the IT model 
(Nevalainen et al. 120031 Kaastra et al. l2003c> . 

To describe the multi-temperature structure in the central 
regions, we used a two temperature model (2T). In Kaastra et 
al. (I2004> and Peterson et al. J2003> we found that in almost all 
cases, at each radius there is negligible emission from gas with 
a temperature less than one-third to half of the ambient temper- 
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(a) Si 



M B7 
A A 1795 
* A 262 



(b) S 



■K- M 87 
A A 1795 
* A 262 



(c) Fe-L 

■iV M B7 
A A 1795 
* A 262 



(d) Fe-K 

■)!r M 87 

A A 1795 

* A 262 



Fig. 1. Comparison between the MOS and pn resuhs for the Si (a), S (b), and Fe (c and d) abundances. We used spectra from 
several regions of M 87 (the brightest cluster at low energies), A 262 (temperature 2 keV), and A 1795 (temperature 4-6 keV) as 
representative samples. 



ature. Therefore, we fixed the temperature separation between 
both components as Tq ~ 2Ti, where Tq and Ti are the two 
temperatures. 

For the inner regions of A 262, 2A 0335+096, A 2052 and 
A 4059, the 2T model improved the fit (TableO. In some cases, 
the Fe abundance changed significantly as compared to the 
single temperature fit (the so-called Fe-bias, see for example 
Buote l2000l and Molendi & Gastaldello 2001 1. For these clus- 
ters, we use the abundances derived from this 2T model for the 
discussion below. For the inner regions in some other clusters 
(e.g. A 496 and A 1795), the 2T model sometimes provides a 
slightly better fit over the IT model. However, in these cases the 
changes in the obtained Fe abundances are small compared to 
the statistical errors. This is because the deviation from isother- 
mality within each radius is relatively small. 

In some cases, the fit is improved further by de-coupling 
To with Ti. However, even in those cases, the change in the 
Fe abundance from the above 2T model is small (~10%) com- 
pared to the errors. 

In all cases except for A 3226, Coma and A 754, the tem- 
perature decreases towards the cluster center (Kaastra et al. 



'2004V In Fig. |5] we present the abundances as derived from 
the acceptable fits (y^ jv < 1.3) of either the IT or 2T mod- 
els. In Tables |3^5] we show average values over three dif- 
ferent spatial scales. We also present average values for cool 
{kTicM < 3 keV), medium-temperature, and hot clusters 
(kTicM > 6 keV) in Table |6l Furthermore, we separately 
show the results for some clusters which exhibit soft excess 
("SB" in TableQ. Fig.|6lshows abundances averaged over (50- 
200) /i^^ kpc for each cluster as a function of the ICM temper- 
ature. Below we summarize the main results of the Fe abun- 
dance, O/Fe, Si/Fe, and S/Fe ratios. 

In all cool and medium-temperature clusters except MKW 
9, the faintest one, we detected a central increase in the 
Fe abundance; the abundance starts to increase at (50- 
100) h^^ kpc in radius from 0.2-0.4 solar to 0.6-0.8 solar to- 
wards the center Except for the Perseus cluster, we found no 
significant Fe increase in the inner regions of the hot clusters 
(See Table|5|l- However, the best-fit value for each cluster sug- 
gests the Fe increase in all cases and the average Fe value over 
the hot clusters also indicates a central enhancement (Table|6|) 
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Table 3. Spectral fit results for the inner regions of some clusters for a single temperature (IT fit) and a two temperature (2T fit) 
model. We present no error when the fit is largely unacceptable. 

cluster W IT fit 2T fit 









Fe/H 


A. 1 ^ 




Fe/H 


A / ^ 






(keV) 


(solar) 




(keV) 


(solar) 




A 262 


0.0-0.5 


1.0 


0.31 


457/227 


1.5 ±0.07 


0.85 ±0.1 


274/226 




0.5-1.0 


1.4 ±0.04 


0.58 ± 0.07 


242/227 


1.6±0.1 


0.89 ±0.1 


214/226 


2A 0335+096 


0.0-0.5 


1.4 


0.32 


400/228 


1.7±0.1 


0.55 ±0.07 


298/227 




0.5-1.0 


1.8 ±0.05 


0.5 ±0.05 


275/227 


2.1 ±0.1 


0.64 ±0.06 


210/227 


A 2052 


0.0-0.5 


1.4 


0.35 


574/228 


2.0 ±0.1 


0.62 ±0.05 


312/227 




0.5-1.0 


2.6 ±0.06 


0.71 ±0.06 


255/228 


2.9 ±0.2 


0.67 ±0.06 


213/227 


A 4059 


0.0-0.5 


2.2 


0.63 


283/226 


2.6 ±0.4 


0.60 ±0.2 


258/223 



^ Radial range (arcmin). 

The best-fit and hotter temperature for the IT fit and for the 2T fit, respectively. 



Table 6. Average abundances obtained from the EPIC analysis for the cool, medium-temperature, hot, and soft excess clusters. 
Same notation as Table 0] 



R 


Fe/H 


0/Fe 


Ne/Fe 


Mg/Fe Si/Fe 


S/Fe 


Ar/Fe 


Ca/Fe 










Cool (fcT < 3 keV) 








#0 


0.57 ±0.04 


0.4 ±0.2 


1.4 ±0.4 


0.6 ±0.3 1.6 ±0.2 


1.4 ±0.3 


1.9 ± 1.5 


3.1 ± 2.1 


#1 


0.42 ± 0.08 


1.9 ± 0.9 


1.0 ±0.8 


- 1.2 ±0.5 


0.7 ±0.6 












medium-T (3 < fcT < 6 keV) 








#0 


0.60 ±0.04 


0.7 ±0.4 


1.8 ±0.5 


0.3 ±0.3 1.6 ±0.3 


1.0 ±0.3 


1.1 ±0.7 


2.9 ± 1.3 


#1 


0.39 ±0.03 


1.2 ±0.5 


1.3 ±0.5 


- 1.4 ±0.4 


0.5 ±0.3 


1.2 ± 1.0 


3.3 ± 1.6 


#2 


0.32 ± 0.08 


4.0 ± 2.5 


1.5 ± 1.2 


- 1.3 ±0.9 


1.7± 1.1 


1.3 ±2.8 












Hot (fcT > 6 keV) 








#0 


0.48 ±0.12 


1.2 ± 1.2 


2.6 ± 1.9 


1.3 ±2.2 2.2 ±1.9 


0.5 ± 1.6 






#1 


0.34 ± 0.05 


2.1 ± 1.7 


3.6 ±2.0 


- 1.9 ±1.7 


1.8 ± 1.8 


3.3 ±4.1 




#2 


0.25 ±0.05 




3.3 ±3.4 




3.1 ±4.0 














Soft excess clusters 








#0 


0.57 ±0.02 


0.6 ±0.2 


1.6 ±0.3 


0.4 ±0.2 1.4 ±0.2 


0.7 ±0.2 


0.8 ±0.4 


2.4 ±0.7 


#1 


0.32 ± 0.02 


0.8 ± 0.3 


0.8 ± 0.4 


- 1.5 ±0.3 


0.4 ±0.3 


1.1 ±0.8 


3.1 ± 1.3 


#2 


0.17 ±0.04 


1.5 ± 1.2 


2.6 ± 1.2 


- 2.2 ±1.0 


1.8 ± 1.1 


1.1 ±3.3 





In the outer regions, (200-500) h^^ kpc in radius for all 
clusters, the Fe abundance is 0.2-0.3, without a significant de- 
pendence on the ICM temperature (Fig.|6j. Exceptions are the 
soft excess clusters, which show relatively smaller Fe abun- 
dances down to 0. 1 solar. These low Fe abundances may not be 
the actual values due to a spectral contamination by the soft ex- 
cess emission, which was not taken into account in the present 
analysis. Further careful analysis is necessary to take into ac- 
count the soft excess emission and derive the accurate Fe abun- 
dance. This is beyond the scope of the present paper. 

The error on the Si abundance determination is larger than 
on iron. Nevertheless, there is also an indication of a central 
increase in the Si abundance. In fact, the Si/Fe ratio is radially 
uniform at an average value of 1-2 solar (Fig.|5j;). Within the 
error bars, there is no radial change of the Si/Fe ratio for each 
cluster and among different clusters (Fig.|6}. The radial distri- 
bution of the S abundance is similar to that of Si; the S/Si ratios 
are consistent with the solar ratio. 

In the cool and medium-temperature clusters, the measured 
O/Fe abundance ratio at the cluster center is 0.4-0.7 relative to 



the solar value. This is consistent with the values obtained from 
the RGS (as shown later). For the hot clusters, the errors on the 
O abundance are too large to allow a measurement of spatial 
variations. Contrary to Fe and Si, the O abundance shows no 
significant central increase. The O/Fe ratio averaged over the 
cool and medium-temperature clusters increases towards large 
radii to a value of a few times solar. 

The Ne, Mg, Ar, and Ca abundances have larger errors. We 
only show the average abundance of these elements in Table |6l 
In general, these abundances relative to Fe are consistent with 
the solar ratio. 

4. The RGS analysis 

4. 1. Analysis method 

The RGS events were filtered using both the dispersion vs. 
cross-dispersion and dispersion vs. pulse height windows. We 
used a spectral extraction region of I'.O in full- width in the 
cross-dispersion direction centered on the X-ray maximum. 
This corresponds to a physical size of 50 h^^ kpc for a mean 
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Fig. 3. An example of the spectral fitting. The spectra of A 262 
between l'-2' in deprojected radius are shown. The MOSl 
and pn spectra are multiplied by factors of 0.25 and 2, respec- 
tively, for clarity. The data are shown with the model as a his- 
togram. Expected line positions for H- and He-like ions are in- 
dicated by vertical lines. The bottom panel shows the fit residu- 
als [(data-model)/model]. In this case, we obtain a temperature 
of 1.9 keV and abundances for O, Si, S, and Fe of 0.5, 1.0, 1.0, 
and 0.7 times solar, respectively. 




0,5 1 2 

Energy (keV) 



Fig. 4. Same as the previous figure, but for the spectra of 
A 1795 between l'-2' in deprojected radius. We obtain a tem- 
perature of 5.3 keV and abundances for O, Si, S, and Fe of 0.20, 
0.7, 0.0, and 0.4 times solar, respectively. 

redshift of 0.06. Generally, most of the source photons in the 
RGS aperture are inside this extraction region. We used a pulse 
height width selection that included 90% of all photons for a 
point source. 

Our targets have a larger spatial extent than the instrumental 
spatial resolution, resulting in a broadening of the line spread 
function (LSF). Here we modeled the source brightness using 
the observed MOS image for each target. Then, we corrected 
for this effect by convolving the cluster brightness profile with 
the RGS response for a point source. Because the dependence 
of the LSF on the off-axis angle within ±3' is smaller than 10- 
30% and the cluster emission is dominated by emission within 
that spatial range, we ignore this dependence of the LSF. 

We have estimated the background spectrum and its vari- 
ability using several blank-sky field observations 400 ks in 
total exposure; Tamura et al. .2003bJ . The count rate variations 



Table 4. Abundances (scaled to the solar ratio) obtained from 
the EPIC analysis for each cluster. Values with errors larger 
than 3 times solar are not presented (indicated by — ). 





Fe/H 


O/Fe 


Si/Fe 


S/Fe 












#0 


0.39 ± 0.06 


0.3 ± 0.7 


1.6 ± 0.7 


2.8 ± 1.3 


#1 


26 ± 09 


1.9 ± 1.9 


1.7 ± 1.2 


1.1 ± 1.7 






A 262 






#0 


67 + n 05 


0.6 ± 0.3 


1.4 ± 0.2 


1.0 ± 0.3 


#1 


0.34 ± 0.04 


1.8 ± 0.9 


1.2 ± 0.4 


0.7 ± 0.5 






A 1837 








n fid + n 1 4 

u.u^ HI U.J- 4: 


n 5 ± 1 fi 


2.8 ± 1.6 


1.0 ± 1.4 


#1 


0.39 ± 0.07 


0.9 ± 1.1 


0.9 ± 0.7 


0.4 ± 0.8 


#2 


27 ± 13 


1.7 ±3.0 


1.0 ± 1.6 


3.7 ± 2.9 






Sersic 159— 


3 




#0 


39 + n 03 


0.6 ±0.4 


1.3 ±0.2 


0.8 ± 0.3 


#1 


0.23 ± 0.03 


0.8 ±0.8 


1.4 ±0.6 


0.5 ± 0.5 


#2 


09 + n 05 




1.5 ± 1.7 


1.5 ± 2.7 






MKW9 






#0 


fi3 + 1 5 


0.1 ±0.5 


1.9 ±0.7 


1.4 ± 0.8 


#1 


0.97 ± 0.45 




0.2 ±0.8 


0.0 ± 0.7 






2A 0335-^096 




#0 


fiO + 04 


1.6 ±0.7 


1.0 ± 0.2 


0.8 ± 0.3 


#1 


49 + fl 07 


1.8 ±1.7 


1.2 ±0.6 


0.0 ± 0.4 


#2 


0.31 ± 0.17 




0.3 ± 1.1 


1.5 ± 1.9 






MKW3s 






#0 


6S + 06 


0.5 ±0.4 


1.4 ±0.3 


0.9 ± 0.4 


#1 


f) oi ^ n no 
vj.OJ- HI xj.yjo 


1.3 ±0.9 


1.5 ±0.6 


n 5 ± n 5 


#2 


0.09 ± 0.06 


0.0 ± 1.9 










A 2052 






#0 


58 + 04 


0.7 ±0.4 


1.3 ±0.3 


0.8 ± 0.3 


#1 


41 + fl 04 


0.6 ±0.5 


1.4 ±0.5 


0.6 ± 0.5 


#2 


0.13 ± 0.06 




0.3 ± 1.1 


3.3 ± 2.6 






A 4059 






#0 


0.72 ±0.12 


0.5 ±0.6 


1.3 ±0.6 


1.7 ±0.8 


#1 


0.44 ± 0.08 


0.7 ±0.7 


1.4 ±0.9 


0.7 ±0.8 


#2 


0.41 ±0.23 


1.2 ±2.5 


3.2 ±2.2 


1.7 ±2.5 






Hyd A 






#0 


0.42 ± 0.05 


0.8 ±0.6 


1.1 ±0.5 


0.4 ±0.5 


#1 


0.30 ±0.05 


0.9 ±0.8 


1.9 ±0.9 


0.7 ±0.8 


#2 


0.20 ±0.11 




0.5 ±1.7 


0.0 ± 1.2 



Integration region (#0, #1, #2 for radii between 0-50, 50-200, 
and 200-500 in h^'^ kpc). 



between the different exposures is small compared to the sta- 
tistical errors: the r.m.s. deviations over the full spectral range 
are less than 30%. Based on this, a 30 % systematic error to the 
background counts was assigned. 

The first order spectra from the RGS 1 and RGS2 were fit- 
ted with the same model simultaneously. We limited the wave- 
length band to 8-28 A, where the estimated background is typ- 
ically less than 10 % of the source counts. Therefore, the sys- 
tematic uncertainty in the background modeUng should not af- 
fect the results significantly. 

In principle, the extracted spectrum contains photons from 
a range of different sky positions. However, in the cases of 
our sample, where the brightness is centrally peaked, the ob- 
served spectrum is dominated by the photons from the cluster 



8 



T. Tamura et al.: The ICM metallicity from XMM-Newton 



Table 5. Abundances obtained from the EPIC analysis for each 
cluster (continued). Same notation as Table |3 



R 


Fe/H 


0/Fe 


Si/Fe 


S/Fe 






A 496 






#0 


0.60 ± 0.05 


0.4 ± 0.3 


1.6 ± 0.3 


1.1 ± 0.4 


#1 


0.37 ±0.07 


1.5 ± 1.2 


1.5 ± 0.8 


0.7 ± 0.8 


#2 


0.41 ± 0.26 




1.7 ± 2.3 


1.8 ± 2.7 






A3112 






#0 


0.70 ± 0.05 


0.9 ± 0.6 


1.4 ± 0.4 


0.4 ± 0.3 


#1 


0.36 ± 0.05 


0.2 ± 0.7 


1.0 ± 0.8 


0.6 ± 0.7 


#2 


0.22 ± 0.13 


0.0 ± 1.6 


1.3 ± 2.3 


2.3 ± 2.7 






A 1795 






#0 


0.51 ± 0.04 


0.7 ±0.4 


1.4 ± 0.4 


0.6 ± 0.4 


#1 


0.32 ± 0.04 


0.5 ± 0.7 


1.9 ± 1.1 


0.0 ± 0.5 


#2 


0.26 ± 0.07 


1.4 ± 1.9 


3.2 ± 2.1 


0.6 ± 1.2 






A 399 






#0 


0.63 ± 0.35 


0.0 ± 2.2 


0.0 ± 1.2 


0.0 ± 1.4 


#1 


0.31 ± 0.12 




0.9 ± 2.2 




#2 


0.28 ±0.13 












A 3266 






#0 


0.46 ± 0.34 








#1 


0.27 ± 0.13 








#2 


0.22 ±0.08 












Perseus 






#0 


0.57 ± 0.03 


1.4 ± 0.3 


1.2 ± 0.2 


0.4 ± 0.4 


#1 


0.41 ± 0.04 


2.3 ± 0.9 


1.2 ± 0.4 


0.2 ± 0.3 






Coma 






#0 


0.53 ± 0.54 








#1 


0.26 ± 0.09 




3.8 ± 2.8 


0.7 ± 2.0 


#2 


0.25 ± 0.10 


2.4 ± 2.8 


1.4± 2.1 


0.0 ± 0.9 






A 754 






#0 


0.32 ±0.34 


0.0 ±2.5 






#1 


0.33 ±0.11 


1.1 ± 2.2 


0.4 ± 1.4 


0.6 ±1.6 


#2 


0.21 ±0.09 












A 1835 






#0 


0.42 ±0.06 


2.0 ± 1.4 


1.1 ±0.9 


0.6 ±1.0 


#1 


0.38 ±0.09 


0.0 ± 1.0 


0.0 ±0.6 


0.0 ±1.1 


#2 


0.29 ±0.11 






1.7 ±2.8 



core. By using a full instrument simulation, SciSim (the official 
XMM-Newton science simulation tool") and the observed MOS 
image, we confirmed that the contamination from the cluster 
outer region (a projected radius from the cluster center larger 
than 3') is less than 10 % for a typical cluster and 20 % at 
most for all clusters. This kind of contamination adds counts 
in the form of a continuum spectrum since lines are smeared 
out due to the large spatial extent of the outer cluster regions. 
This could make the obtained equivalent line width and hence 
the metal abundance relative to H appeared lower than the ac- 
tual value by less than 10-20 %. The possible change of abun- 
dances relative to Fe (e.g. O/Fe) should be much smaller Since 
the thermal structure is mainly determined by the Fe-L line ra- 
tios, the effect of the contamination on the thermal structure 
should also be small. Therefore, we ignored this flux mixing 
and assumed that the values obtained from our spectral fits are 
emission-weighted values for the cluster core. 



< 
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Fig. 6. Abundance vs temperature for each cluster averaged 
over (50-200) h^^ kpc. Errors on O/Fe and S/Fe are not shown 
for clarity (These are presented in Tables |4j|5ji- The data for 
MKW 9 are not shown for clarity, because this cluster has 
larger errors than the other clusters. 



4.2. Single-temperature model fitting 

We fitted a set of spectra for each cluster with the single tem- 
perature (IT) model. With this case, the free parameters are the 
normalization, kT, the Fe abundance, and the Fe/O ratio. An 
exception is NGC 533, where the RGS spectrum is dominated 
by line emission. In this case, it is relatively difficult to deter- 
mine the continuum level and hence the line equivalent widths. 
Therefore for NGC 533 we fixed the Fe abundance to the value 
obtained from the EPIC analysis above. 

We obtained acceptable or almost acceptable fits (y^ jv < 
1.2) with the IT model for Sersic 159-3, A 4059, Hyd-A, 
A 496, A 3112, A 1795 and A 1835. In contrast, NGC 533 
(x^/v = 1.55), A 262 (1.49), A 2052 (1.32), and 2A 0335 
(1.48) exhibit clear deviations from the IT model. 

4.3. Multi-temperature model fitting 

In some clusters, there is evidence for multi-temperature emis- 
sion (Peterson et al.,2003 l. If there is a significant contribution 
from the cooler temperature components in the Une emission, 
the abundances derived from the IT model may differ from the 
true abundances. 

To avoid this error, we used a multi-temperature model 
which consists of at most five temperature CIE components 
with fixed temperature separation, (Tq = 2Ti = 4T2 = 8T3 = 
I6T4), but free normalization (EMi) for each component (the 
5T model); a similar approach was also taken by Peterson et 
al. (2003 1. In this case, the Ne/Fe and Mg/Fe ratios are also 
allowed to be free along with the Fe/H and O/Fe ratios. All 
temperature components were assumed to have the same abun- 
dances. 
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Table 7. RGS fit results for a multi-temperature model. The mean abundances for the sample are also shown in the last row. 



target 




kTS 


O/Fe 


Ne/Fe 


Mg/Fe 


Fe/H 




NGC 533 


1, 


.4 ±0.5 


0.54 ±0.15 


1.3 ±0.7 


1.7 ± 1.0 


0.50 


264/226 


A 262 


2, 


.4 ±0.7 


0.64 ±0.12 


0.6 ±0.5 


0.1 ±0.5 


0.58 ±0.11 


229/231 


Sersic 159-3 


2, 


.3 ±0.2 


0.67 ±0.15 


0.8 ±0.4 


0.6 ±0.5 


0.38 ±0.03 


279/263 


A 2052 


2, 


.8 ±0.3 


0.66 ±0.10 


0.8 ±0.3 


0.9 ±0.6 


0.68 ±0.08 


284/263 


2A 0335 


1, 


.6 ±0.1 


1.00 ±0.23 


0.4 ±0.5 


1.0 ±0.6 


0.29 ±0.04 


282/188 


A 4059 


3 


.5 ±0.7 


0.68 ±0.16 


1.4 ±0.7 


0.6 ±0.8 


0.58 ±0.10 


215/181 


Hydra A 


2, 


.9 ±0.3 


0.70 ±0.18 


0.5 ±0.6 


0.0 ±0.3 


0.37 ±0.06 


230/228 


A 496 


3, 


.2 ±0.6 


0.66 ±0.15 


1.2 ±0.6 


1.1 ±0.9 


0.55 ±0.14 


236/229 


A3112 


4, 


.2 ±0.8 


0.63 ±0.18 


0.5 ±0.6 


0.0 ±0.4 


0.69 ±0.14 


189/225 


A 1795 


6, 


.4± 1.0 


0.78 ±0.14 


1.5 ±0.4 


1.4 ±0.8 


0.60 ±0.13 


284/261 


A 1835 


4, 


.6± 1.0 


0.35 ± 0.24 


0.9 ±0.6 


0.0 ±0.7 


0.49 ±0.12 


295/284 


mean 






0.63 ±0.05 


1.0 ±0.2 


0.66 ±0.2 


0.54 ±0.03 





^ Temperature of the hottest component in keV. 



The results are shown in Table0 In Sersic 159—3, A 4059, 
Hyd-A, A 3112 and A 1835 where the IT model already pro- 
vided an acceptable fit, this model does not improve the fit. The 
derived abundances are consistent with those of the IT fit, sim- 
ply because the same temperature component dominates the 
emission in both the IT and 5T models. 

In 2A 0335H-096, we could not obtain an acceptable fit even 
with the 5T model. We note that the EPIC spectra of the inner 
region of this cluster also show a deviation from the 2T model 
(Table|3ll. The Galactic absorption towards this cluster is by far 
larger than that of other clusters. Furthermore, the X-ray image 
shows a complex structure (Sarazin et al. 1992 1. Therefore, we 
assumed that our spectral and spatial modeling is oversimpli- 
fied for this cluster. Note that we exclude the EPIC results for 
the inner region of this cluster in our discussion. 

In other clusters, the 5T model improves the fit signifi- 
cantly. The best example is A 262, where the model reproduced 
the observed spectra very well. By using this model, we ob- 
tained acceptable fits jv < 1.2) for all clusters. While there 
are significant changes in the derived Fe abundance in some 
cases, the change in the O/Fe ratio is at most 30 % or 0.2 times 
solar, compared to the IT fitting. 

Peterson et al. ( I2003t analyzed a similar RGS data set, fo- 
cusing on the temperature structure. Our model for the temper- 
ature structure is basically consistent with their results. 

In Fig. 13 we compare the O and Fe abundance as obtained 
by EPIC and RGS. In 2A 0335H-096, the Fe/H ratio derived 
from the RGS data is two times smaller than that derived from 
the EPIC data, while the O/Fe ratio from the two instruments 
is consistent. We presume that this discrepancy is due to the 
nature of the central spectra which is more complex than our 
modeling. Although the errors in the oxygen abundance deter- 
mination from EPIC is large, apart from 2A 0335H-096, the two 
instruments provide almost consistent results: a formal fit of 
the Fe abundance from RGS on the Fe abundance from EPIC 
to the 1:1 relation yields x^/i^ = 7.68/10. 



5. Summary and Discussion 

5. 1. Present results 

Using spatially-resolved EPIC spectra, we have resolved the 
temperature and abundance distribution on a sub-arcmin scale 
in a sample of 19 clusters. Furthermore, using high spectral res- 
olution RGS spectra of 1 1 cluster cores, we have resolved the 
soft X-ray spectrum dominated by the Fe-L and the O VIII Lya 
lines. This provides not only tight constrains on the thermal 
structure but also a robust measurement of the O abundance. 
Below we summarize our main results. 

1 . We found no significant variation in the Fe abundance and 
the Si, S, and O ratios to the Fe among the systems (see 
Fig.EJ. 

2. The O/H and O/Fe ratio in the cluster cores, (20- 
60) h'^ kpc in radius, are 0.34±0.03 and 0.63±0.05, re- 
spectively. The rm.s. deviation for individual clusters as 
compared to these average values is comparable to the mea- 
surement errors on these ratios. 

3. In all clusters with a temperature less than 6 keV, except for 
the faint cluster MKW 9, we detected a central increase in 
the Fe abundance; the abundance starts to increase at (50- 
100) h^^ kpc in radius from 0.2-0.4 solar value to 0.6- 
0.8 solar value towards the center. 

4. The Si/Fe and S/Fe ratios in cool and medium tempera- 
ture clusters are radially uniform within (200-500) kpc 
with mean values of 1 .4±0.2 and 1 . 1 ±0.3, respectively, rel- 
ative to the solar ratio. 

5. Contrary to the Fe, Si, and S abundances, the O abundance 
shows no spatial variations. The errors on the O abundance 
for each cluster are large. However, when we combine the 
results from several clusters and assume no correlation of 
the systematic errors of the measurements between clusters, 
we detected a significant radial variation in the O/Fe ratio; 
the mean O/Fe ratios are 0.5to'i l-7±04' and A.ltH for 
(0-50) h-^ kpc,(50-200) h-^ kpc, and (200-500) h-^ kpc, 
respectively. 
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5.2. Similarity of the abundance among clusters 

ASCA and BeppoSAX observations have established the Fe in- 
crease towards the cluster center in relatively cool clusters with 
cD galaxies (e.g. Fukazawa et al. 120001 De Grandi and Molendi 
1200 1 1. Our measurement not only confirms these results in the 
cool and medium temperature clusters but also suggests a simi- 
lar enhancement in hot clusters without a prominent cD galaxy 
(e.g. Coma), probably because of the improved spatial resolu- 
tion. 

Beside the Fe abundance, we found no variation in the Si 
and other element ratios to Fe among the clusters. This im- 
plies a homogeneous production of the metals in our sample. 
However, the sample is biased towards medium-temperature 
X-ray luminous clusters. Therefore we make no conclusion 
about the universal distribution of the relative abundances. In 
fact Fukazawa et al. (Jl998)l claimed an increase in Si/Fe from 
cool to hot systems using 40 clusters. 

5.3. The O abundance at the cluster core 

The theoretical models for supernova (SN) types la and II pre- 
dict the O/Fe ratio to be < 0.05 and 1.5 — 4 times solar val- 
ues, respectively. On the other hand, the observed O/Fe ratio at 
the cluster cores are between these two predicted values. This 
suggests that the ICM metals at the cluster cores could be pro- 
duced by a mix of both types of SNe. A value of the O/Fe abun- 
dance ratio in the ICM close to solar values supports the idea 
that the metals in the ICM and the Galaxy were produced in 
a similar way. Assuming the O and Fe yields of SN la and II 
from Tsujimoto et al. ( 1995 ), the observed O/Fe ratio (^ 0.65) 
can be reproduced by the relative frequencies of the two SNe, 
Nin/Nii, of ^ 0.6. In this model, the mass fractions of the O 
and Fe originated from SN la are 0.05 and ^ 0.8, respectively. 

The observed O abundance provides an accurate way to de- 
termine the total mass of O in the ICM, M}^^. We estimate 
M}p^ within 50 h-^ kpc to range from 10** (NGC 533) to 
2 X 10^ (A 1795) in units of /I'^.s ^j. q^j. rq^ sample. 
We assume that all oxygen was produced in SNe II. Then, com- 
bining M}^^ with a SN II model O yield (- 2 Mq, Tsujimoto 
et al.[T995l. gives a total number of SNe II of 10^ - 2.5 x 10^ 
for h = 0.7. Further assuming that the duration of SNe II ac- 
tivity is 10*^ yr ( a typical life time of a 20 Mq star), we can 
estimate the SN II rate to be 10-200 SN He yr"^ This rate is 
at least one or two orders of magnitude larger than that of a 
typical starburst galaxy. In other words, the observed O mass 
requires a large amount of massive star formation, producing 
the corresponding number of SNe II. 

The observed O/H ratio in cluster cores (20-60 kpc in ra- 
dius) is about half that observed in the central regions (5- 
15 kpc) of giant elliptical galaxies, like in M87 (Sakelliou et al. 
I2OO2I Gastaldello and Molendi 2002 Finoguenov et al. l2002l l 
and NGC 4636 (Xu et al. 2002 1. The overabundance of oxy- 
gen in galaxy centers suggests that in those regions oxygen has 
been produced via massive star formation and that metals have 
not been completely mixed with the ambient gas (Gastaldello 
and Molendi .2002j . 



5.4. The total Fe mass in clusters 




10 100 

R (ir' kpc) 



Fig. 7. The Fe number density in the ICM. Same notation as 
Fig. 13 
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Fig. 8. Same as the previous figure, but for the integrated Fe 
mass in the ICM. 



Table 8. The iron mass to light ratio within 250 h ^ kpc. 



cluster 


Mfc 




Mfc/Lb 




(10» Mo) 




(IO-^Mq/Lo) 


MKW9 


0.5 


13 


0.3 ±0.15 


MKW3s 


1.5 


13 


0.9 ±0.2 


A 496 


2.8 


14 


1.5 ±0.5 


A 1795 


6.0 


27 


1.0 ±0.2 


A 754 


3.8 


29 


0.7 ±0.3 



Number of galaxies within a radius of 250 h ^ kpc. 
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We also calculated the total Fe density and mass in the 
ICM as a function of radius, based on our measurements of 
the emission measure and the Fe abundances (Figs. and 
To compare this Fe mass and the predicted Fe production rate 
in galaxies, we estimate the iron mass to light ratio (IMLR). 
We use the galaxy number within a radius of 250 h^^ kpc, 
-/Vo.5 (Bahcall 1981 and reference therein), and estimate the to- 
tal optical luminosity of member galaxies based on the rela- 
tion Lb = 4.5 X 10^N^:^^h-^ Lq (Edge and Stewart 1991 1. 
Table [S] shows the results for representative clusters, where 
-/Vo.5 is available. 

The IMLR of MKW 9 is relatively smaller than for other 
clusters. The potential of MKW 9 might not be as deep as in 
other systems and a large part of the metal-enriched gas has 
been driven out already from the potential well. An IMLR of 
order 10^^ in other systems is consistent with previous mea- 
surements (e.g. Tsuru lTWrb and at least a few times larger than 
the predicted value based upon the current rate of SN la in ellip- 
tical galaxies, as argued previously (e.g., Tsuru .l991t Renzini 
et al. ll993> . We presume that a large fraction of Fe has been 
produced also in the early phase of the galaxy formation with 
a high production rate. Relatively low metal abundances ob- 
served in elliptical galaxies (sub-solar at most; e.g. Matsushita 
et al. EOOOt also support this idea. 

5.5. Spatial distribution of t lie O, Si, S, and Fe 
abundances 

The most simple model to explain the results (3) and (4) is that 
Si, S, and Fe have a common origin for a large part, while O has 
a different origin. This contrast between Si-S-Fe and O has re- 
ported separately in A 496 (Tamura et al. 2001b). Qualitatively 
this trend is consistent with the current idea that the contribu- 
tion from SN la relative to that from SN II, A^ia/^ii, increases 
towards the cluster center and the O/Fe ratio is more sensitive 
to Nin/Nii than the Si/Fe ratio. We do not attempt a further 
detailed comparison between the observation and models, be- 
cause of large uncertainties in both the data (e.g. the O/Fe ratio 
in the outer regions and for hot clusters) and the theory (e.g. 
Gibson et al. .l997i ). 

Acknowledgements. This work is based on observations obtained with 
XMM-Newton, an ESA science mission with instruments and contri- 
butions direcfly funded by ESA Member States and the USA (NASA). 
The Space Research Organisation of the Netherlands (SRON) is 
supported financially by NWO, the Netherlands Organisation for 
Scientific Research. We thank the referee for useful comments. 



Finoguenov, A., David, L.R, Ponman, T.J., 2000, ApJ, 544, 188 
Finoguenov, A., Matsushita, K., Bohringer, H., Ikebe, Y, & Arnaud, 

M. 2002, A&A, 381, 21 
Fukazawa, Y., Makishima, K., Tamura, T., et al. 1998, PASJ, 50, 187 
Fukazawa, Y, Makishima, K., Tamura, T, Nakazawa, K., Ezawa, H., 

Ikebe, Y, Kikuchi, K., & Ohashi, T. 2000, MNRAS, 313, 21 
Gasteldello, R & Molendi, S. 2002, ApJ, 572, 160 
Gibson, B. K., Loewenstein, M., & Mushotzky, R. F. 1997, MNRAS, 

290, 623 

De Grandi, S. and Molendi, S., 2001, ApJ, 551, 153 
den Herder, J.W., Brinkman, A.C., Kahn, S.M., et al. 2001, Ac&A, 365, 
L7 

Jansen, R, Lumb, D., Altieri, B., et al. 2001, A&A, 365, LI 
Kaastra, J.S. 1992, An X-Ray Spectral Code for Optically Thin 

Plasmas (Internal SRON-Leiden Report, updated version 2.0) 
Kaastra, J.S., Ferrigno, C, Tamura, T, et al. 2001, A&A, 365, L99 
Kaastra, J.S., Lieu, R., Tamura, T, Paerels, F.B.S., den Herder, J.W. 

2003a, A&A, 397, 445 
Kaastra, J.S., Tamura, T, Peterson, J.R., et al., 2004, A&A, 413, 415 
Kaastra, J.S., Lieu, R., Tamura, T, Paerels, F.B.S., den Herder, J.W. 

2003c, Adv. Sp. Res., in press. 
Liedahl, D.A., Osterheld, A.L., & Goldstein, W.H. 1995, ApJL, 438, 

115 

Matsushita, K., Ohashi, T, & Makishima, K. 2000, PASJ, 52, 685 
Matsushita, K., Finoguenov, A., & Bohringer, H. 2003, A&A, 401, 
443 

Mewe, R., Gronenschild, E.H.B.M., and van den Oord, G.H.J. 1985, 

A&AS, 62, 197 
Molendi, S., & Gastaldello, R, 2001, A&A, 375, L14 
Nevalainen, J., Lieu, R., Bonamente, M., & Lumb, D., 2003, ApJ, 584, 

716 

Peterson, J.R., Paerels, F.B.S., Kaastra, J.S., et al., 2001, A&A, 365, 
L104 

Peterson, J.R., Kahn, S.M., Paerels, RB.S., et al. 2003, ApJ, 590, 207 
Renzini, A., Ciotti, L., D'Ercole, A. & Pellegrini, S., ApJ, 1993, 419 
Sarazin, C. L., O'Connell, R. W., & McNamaia, B. R. 1992, ApJL, 
397, L31 

Sakelliou, I., Peterson, J.R., Tamura, T, et al. 2002, A&A, 391, 903 
Struder, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18 
Tamura, T, Kaastra, J.S., Peterson, J.R., et al., 2001a, A&A, 365, L93 
Tamura, T, Bleeker, J.A.M., Kaastra, J.S., et al. 2001b, A&A 379, 107 
Tamura, T, den Herder, J.W., & Gonzalez-Riestra, R., 

2003, XMM calibration note, XMM-SOC-CAL-TN-0034 

I http://xmm.vilspa.esa.es/ 1 
Tsujimoto, T, Nomoto, K., Yoshii, Y, Hashimoto, M., Yanagida, S., 

& Thielemann, R-K. 1995, MNRAS, 277, 945 
Tsuru, T. 1991, PhD thesis. University of Tokyo 
Turner, M.J., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27 
Xu, H. et al. 2002, ApJ, 579, 600 



References 

Anders, E., Grevesse, N., 1989, Geochim. Cosmochim. Acta 53, 197 
Arnaud, M., Rothenflug, R., Boulade, O., Vigroux, L., & Vangioni- 

Flam, E. 1992, A&A, 254, 49 
Bahcall, N. A. 1981, ApJ, 247, 787 
Buote, D.A., 2000, MNRAS, 311, 176 

David, L. P., Nulsen, P. E. J., McNamara, B. R., Forman, W., Jones, 

C., Ponman, T, Robertson, B., & Wise, M. 2001, ApJ, 557, 546 
Dickey, J.M., & Lockman, RJ., 1990, ARA&A 28, 215 
Edge, A.C., & Stewart, G.C. 1991, MNRAS, 252, 428 



12 



T. Tamura et al: The ICM metallicity from XMM-Newton 

o 



(a) Fe/H 



X 



in 
d 



X 



d 




□ 



X 



10 100 
R (h^^ kpc) 



(b) 0/Fe 




10 100 
R (h^^ kpc) 



(c) Si/Fe 



C) 




(d) S/Fe 



C) 



o 



C) 

^- 

XX 



X 




o 



10 



100 



o 



10 



100 



R (h ^ kpc) 



R (h ^ kpc) 



Fig. 5. Radial profiles of the Fe/H, 0/Fe, Si/Fe and S/Fe ratios in solar units, derived from the EPIC spectra. The best-fit values 
for cool, medium-temperature, hot, and soft excess clusters are shown with open-circles, filled-circles, open-square, and "x", 
respectively. Typical error bars are shown for A 262, A 496, and Perseus. Mean values and errors for the cool and medium- 
temperature clusters (T< 6 keV), except for the soft excess clusters, are shown with diamonds. In panel (b), the O/Fe ratio 
obtained from the RGS spectra are also shown with fiUed-stars. 



